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ABSTRACT 


The  multiobject  tracking  problem  is  one 
tent  ion  iri  recent  years.  The  primary  difficulty  i 
dating  measurements  among  ^ several  tiumfrorfclperhc 
In  the  present  work,  a  recently  proposed  algorithm 
I'or  application  in  distributed  processing  passive 
by  distributing  the  measurement- target  association 
al  sites  in  the  network  and  then  combining  each  si 
tions  to  form  the  desired  target  tracks.  Compute 
demonstrate  the  capabilities  of  the  procedure. 
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CHAPTER  1 
INTRODUCTION 


The  need  for  ever  more  sophisticated  multi-target  algorithms  has  in¬ 
creased  greatly  in  recent  years.  In  military  applications,  especially,  the 
need  for  these  techniques  is  evident.  Recently,  there  has  been  increased 
interest  in  distributed  tracking  systems,  in  which  the  tracking  problem  is 
broken  down  into  several  smaller  problems  and  distributed  to  trackers  at 
several  different  sites.  Such  problems,  of  course,  present  more  challenges 
than  the  more  traditional  central  processing  tracking  systems. 

The  basic  problem  in  multi-target  tracking  is  data  association.  If 
there  are  several  (perhaps  an  unknown  number  of)  targets  in  an  area,  meas¬ 
urements  are  received  from  each  of  them.  There  exists  an  uncertainty  in 
the  true  origin  of  each  measurement.  Essentially,  any  multi-target  track¬ 
ing  procedure  should  partition  the  measurements  into  sets  associated  with 
each  target  before  the  tracking  of  that  target  can  be  done. 

Several  important  papers  that  have  appeared  in  recent  years  come  to 
grips  with  thi3  problem  in  a  variety  of  ways.  Bar-Shaloin  L3J  gives  an  ex¬ 
cellent  survey  of  recent  work  in  the  field.  The  most  important  paper  for 
this  work  is  the  one  by  Donald  Reid  [S3.  Reid  formulates  hypotheses  of  the 
origin  of  measurements  sequentially,  and  organizes  them  in  the  form  of  a 
hypothesis  tree.  For  each  hypothesis  constructed,  a  set  of  Kalman  filters 
is  constructed  to  track  the  targets  implied  by  that  hypothesis.  The 
strength  of  the  hypothesis  is  then  evaluated  by  comparing  the  actual  meas- 
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uremeuts  with  predicted  values  obtained  from  the  state  estimates  of  the 
filters.  A  specific  formula  for  this  is  developed  by  Reid,  and  incor¬ 
porates  false  alarms  and  new  targets. 

Reid  assumes  in  his  work  a  central  processing  system  -  that  is,  a  sys¬ 
tem  in  which  all  tracking  computations  are  done  at  a  single  site.  The  goal 
of  this  work  is  to  apply  Reid's  ideas  to  a  distributed  processing  system. 
As  we  shall  see,  the  data  association  and  target  tracking  can  be  separated 
quite  naturally  to  fit  in  the  mold  of  such  a  system.  Our  primary  example 
will  be  a  two  node  passive  system  that  receives  bearings  to  targets  as 
measurements.  The  information  signal  is  assumed  to  be  acoustic  in  nature, 
which  in  combination  with  a  desire  to  track  fast  moving  targets  improves 
the  problems  of  propagation  delay. 

Chapter  2  develops  the  theoretical  structure  of  the  proposed  tracking 
algorithm  in  great  detail.  In  Chapter  j,  we  discuss  the  concept  of  a  Dis¬ 
tributed  Sensor  Network  (DSN),  which  motivates  the  choice  of  our  primary 
example.  Chapters  **  and  5  develop  the  mathematical  formulations  necessary 
for  the  implementation  of  the  tracking  algorithm.  Finally,  in  Chapter  6, 
we  present  some  demonstrations  of  the  performance  of  the  algorithm. 


CHAPfER  2 


THE  ALGORITHM 

In  this  chapter,  we  will  model  the  solution  of  a  generalized  tracking 
problem,  and  see  how  previous  work  in  multi-target  tracking  fits  into  this 
model . 

2.1  SOME  DEFINITIONS 

Consider  a  generalized  tracking  system  T.  The  task  of  T  is  to  deter¬ 
mine  the  existence  of  and  track  certain  objects  located  in  a  specified  en¬ 
vironment  E.  The  objects  are  called  targets .  The  system  tracks  a  target 
if  it  can  estimate  its  position  and  velocity  at  any  time  while  it  is  locat¬ 
ed  in  E.  A  tracking  system  may  also  -determine  acceleration  and  other 
higher  positional  derivatives;  however,  for  our  purposes,  position  and 
velocity  will  be  sufficient.  This  information  is  collectively  called  the 
target  state. 

To  perform  its  task,  the  tracking  system  must  obtain  information  from 
the  environment.  Usually,  this  information  is  in  a  form  termed  signals. 
For  example,  a  radar  system  might  use  electromagnetic  sinusoids,  whereas  in 
sonar,  the  signals  are  acoustic.  Any  place  within  T  that  receives  signals 
is  called  u  sensor.  If  the  received  signals  originated  in  T  and  were  re¬ 
flected  back,  T  is  called  an  active  system.  If  the  signals  originated  in 
the  environment,  T  is  a  passive  system.  It  may  happen  that  the  system  re¬ 
ceives  both  types  of  signals,  in  which  case  'i'  is  referred  to  a3  a  mixed 
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system.  Active  systems  generally  obtain  more  information  from  the  environ¬ 
ment  than  passive  systems;  however,  they  are  also  more  susceptible  to 
detection  by  the  targets  themselves.  A  tracking  system  can  also  be  classi¬ 
fied  by  the  manner  in  which  it  collects  signals.  Our  attention  will  focus 
on  time-sampling  systems,  which  sample  signals  from  all  directions  at 
specified  instants  of  time. 

Usually,  the  received  signals  in  their  original  form  are  not  very  use¬ 
ful  for  tracking.  A  signal  processor  in  T  extracts  data  from  the  signals 
which  can  be  used  in  tracking.  In  this  work,  data  obtained  in  this  manner 
from  signals  are  termed  measurements.  Measurements  may  be  scalar  or  vector 
in  nature.  For  example,  a  measurement  from  a  radar  system  might  consist  of 
three  position  coordinate?1.  A  measurement  from  a  passive  system,  on  the 
other  hand,  might  simply  be  a  bearing.  In  time-sampling  systems,  a  scan  is 
j  defined  as  a  set  of  measurements  obtained  at  the  same  time-sampling  in- 

!  stant.  In  many  systems,  a  signal  processer  requires  signals  from  several 

different  sensors  to  produce  a  single  measurement.  Normally,  this  set  of 
sensors  is  fixed  and  is  referred  to  collectively  as  a  node.  Since  our  at- 

! 

tention  is  focused  on  the  tracking  problem,  we  will  not  go  into  detail  the 
highly  nontrivial  problem  of  signal  processing.  Hence,  all  of  our  work 
will  refer  to  measurements  and  nodes  at  the  lowest  level. 

/ 

j 

!  2.2  -  THE  TRACKING  PROBLEM 


In  general,  the  tracking  system  must  be  able  to  track  several  targets 
simultaneously.  In  multi-target  tracking,  there  are  really  two  problems 
that  must  be  solved.  They  are  the  data  association  problem  and  the  track 
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association  problem .'  Logically,  before  a  tracker  can  estimate  the  state  of 
a  particular  target,  it  must  know  which  measurements  to  use.  Given  a 
time-sequenced  set  of  scans,  the  tracker  must  partition  the  set  of  all 
measurements  into  a  set  of  mutually  exclusive  subsets.  One  subset  may 
correspond  to  an  assignment  to  no  target  at  all.  The  measurements  in  this 
set  are  called  false  alarms.  The  other  subsets  will  correspond  to  hy¬ 
pothesized  targets.  The  process  of  forming  these  subsets  is  data 
association,  and  each  subset  of  measurements  is  called  a  data  track. 

There  will  be  many  different  ways  to  perform  data  association  on  a  given 
set.  Each  such  configuration  is  called  a  data  hypothesis. 

We  will  asssume  in  our  work  that  scans  are  to  be  processed  sequential¬ 
ly.  This  is  reasonable,  since  in  a  real-time  system,  it  is  usucilly  desir¬ 
able  to  incorporate  new  information  as  soon  as  possible.  Therefore,  at  any 
given  time,  the  number  of  possible  data  hypotheses  depends  on  the  number  of 
a  priori  hypotheses  and  the  number  of  measurements  in  the  present  scan. 
The  sequential  formation  of  hypotheses  can  be  organized  conceptually 
through  the  use  of  trees. 

The  structure  of  hypothesis  trees  is  best  illustrated  through  a 
specific  example.  Assume  a  priori  there  are  no  data  tracks.  At  time  tsO, 
the  tracker  receives  two  measurements.  We  make  the  simplifying  assumption 
that  a  particular  target  cannot  be  the  source  of  more  than  one  measurement 
in  a  given  scan,  although  the  more  general  case  can  still  be  displayed  in 
tree  form.  Thus,  each  measurement  could  be  from  a  different  legitimate 
target  or  a  false  alarm.  This  results  in  four  different  data  hypotheses, 
as  shown  in  Fig.  2-1. 
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ml  and  m2  are  the  two  measurements.  As  can  be  seen,  each  corresponds 
to  o  different  level  of  the  tree.  The  locations  within  the  tree 
corresponding  to  the  assignment  of  measurements  to  a  data  track  are  called 
nodes.  (The  use  of  this  terms  is  unfortunate;  however,  it  will  always  be 
clear  from  oontext  whether  we  are  talking  about  groups  of  sensors  or  tree 
structure.)  A  sequence  of  connected  nodes  is  called  a  branch  of  the  tree. 
Each  branch  of  the  tree  corresponds  to  a  different  data  hypothesis,  as 
shown.  Thus,  hypothesis  1  assigns  ml  to  data  track  1  and  m2  to  data  track 
2,  while  hypothesis  3  assigns  ml  as  a  false  alarm  (represented  by  0)  and  m2 
to  data  track  2.  Ue  have  made  the  convention  of  identifying  a  data  track 
by  the  measurement  number  of  its  first  measurement. 

Assume  now  that  we  have  a  second  scan  with  two  measurements  (m3  and 
m4).  The  number  of  assignments  of  m3  and  m4  will  depend  on  which  prior  hy¬ 
pothesis  is  assumed.  For  instance,  if  hypothesis  1  is  assumed,  the  possi¬ 
ble  assignments  fo-  each  are  data  track  1,  data  track  2,  a  new  data  track, 
or  a  false  alarm,  witn  the  condition  that  both  cannot  be  assigned  to  the 
same  legitimate  data  track.  If  hypothesis  4  is  assumed,  the  only  possible 
assignments  are  a  new  data  track  or  a  false  alarm.  If  we  expand  each 
branch  .  in  this  manner,  we  obtain  Fig.  2-2,  which  represents  the  total 
number  of  ways  the  measurements  in  the  two  scans  can  be  associated,  given 
our  assumptions.  It  is  easy  to  see  that  the  tree  expands  at  an  exponential 
rate . 


The  above  procedure  for  the  construction  of  the  hypothesis  tree  is 
known  as  a  measurement-oriented  approach,  because  every  possible  data  track 


is  listed  for  each  measurement 


A  target-oriented  approach  would  list 


every  possible  measurement  for  each  data  track.  However,  in  taking  the 
latter  approach,  it  is  conceptually  difficult  to  decide  when  a  new  data 
track  should  be  hypothesized,  whereas  it  .ie  former  case,  new  datu  tracks 
appear  naturally  as  a  part  of  tree  expansion.  For  this  reason,  we  will  use 
the  measurement-oriented  approach. 

If  there  are  mul  iple  nodes  in  the  tracking  system,  we  oan  perform  the 
data  association  in  two  ways.  A  central  processing  system  would  construct 
one  hypothesis  tree,  incorporating  all  measurements  from  every  node  in  the 
system.  A  distributed  processing  system  would  construct  a  tree  for  each 
node  or  a  subset  of  nodes.  We  can  think  of  nesting  the  various  trees  in  a 
distributed  processing  system  in  one  another  to  produce  an  overall  data  as¬ 
sociation  hypothesis  tree  that  is  equivalent  to  the  tree  constructed  in  a 
central  processing  system,  which  contains  all  possible  data  tracks  that  can 
be  formed  given  the  received  measurements.  With  a  distributed  structure, 
we  encounter  the  track  association  problem. 

The  track  association  problem  can  also  be  modeled  as  a  hypothesis 
tree.  Whereas  data  hypotheses  associate  measurements,  track  hypotheses  as¬ 
sociate  the  data  tracks  of  different  nodes.  Given  a  track  hypothesis  and 
the  data  tracks  it  is  conditioned  on,  we  finally  can  combine  the  data  to 
form  sets  of  estimated  target  states,  or  target  tracks .  (This  assumes,  of 
course,  that  there  are  a  sufficient  number  of  data  tracks  to  make  sucli  on 
estimation  possible;  i.e.  the  system  must  be  observable) .  The  production 
of  these  target  tracks  is  the  desired  goal  of  the  tracking  system.  Since 
the  track  association  trees  are  each  conditioned  on  data  hypotheses  from 
the  various  nodes  of  the  system,  we  can  nest  the  former  into  the  overall 


data  association  tree  of  the  tracking  system.  The  resulting  structure  or¬ 
ganizes,  in  a  systematic  way,  ail  possible  solutions  of  the  tracking  prob¬ 
lem. 

The  prospect  of  wading  through  such  an  imposing  structure  is  horrify¬ 
ing,  at  best.  In  the  next  sections  we  will  discuss  ways  of  making  the 
solution  more  tractable.  These  methods,  of  course,  produce  suboptimal 
results,  in  the  sense  that  the  correct  solution  may  accidentally  be  dis¬ 
carded.  However,  they  are  quite  necessary  for  any  oractical  implementa¬ 
tion  . 

£•1  “  DATA  HYPOTHESIS  REDUCTION  -  TREE  PRUNING 

The  use  of  trees  as  a  problem-solving  technique  is  well  knowr.  in  ar¬ 
tificial  intelligence.  They  arc  used  to  systematically  model  the  step-by- 
step  solution  of  very  general  problems.  As  we  have  seen,  they  fit  natural¬ 
ly  into  the  tracking  problem.  However,  in  most  applications,  there  is  a 
well-defined  "goal  state",  which  will  be  reached  eventually  by  one  of  the 
branches  of  the  tree.  In  the  tracking  problem,  there  is  no  "goal  state"; 
the  trees  are  open-ended.  Since  the  trees  expand  exponentially,  we  must 
impose  constraints  on  tree  growth  in  order  to  keep  the  problem  manageable. 
Methods  U3ed  for  this  purpose  arc  called  tree  pruning  techniques. 

There  are  two  ways  to  limit  tree  expansion.  The  breadth  of  the  tree 
can  be  limited  by  retaining  no  more  than  a  specified  maximum  number  of 
branches  each  time  the  tree  expands.  This  "pruning"  of  branches  of  the 
tree  may  result  in  one  branch  of  an  older  scan  being  singled  out.  This  is 


illustrated  in  Fig.  2-3.  In  this  case,  the  scan  is  said  to  be  identified. 


Scans  will  not  always  be  identified  by  breadth  constraints.  Therefore,  we 
must  also  limit  the  depth  of  a  tree.  When  the  number  of  scans  in  a  tree 
attains  a  certain  maximum  allowable  number,  one  branch  in  the  oldest  scan 
is  singled  out  and  the  others  are  pruned.  This  scan  is  thus  forcibly  iden¬ 
tified.  In  limiting  the  depth  of  a  tree,  redundant  hypotheses  may  appear 
in  the  tree,  depending  on  the  mechanism  used.  A  set  of  redundant  hy¬ 
potheses  in  a  depth  limited  tree  assign  measurements  in  the  same  manner. 
They  may  assign  measurements  to  exactly  the  same  targets  (in  which  case 
they  are  termed  identical)  or  there  may  exist  a  1-1  mapping  of  targets 
between  the  hypotheses.  A  set  of  redundant  hypotheses  may  be  combined 
under  the  assumption  that  any  past  differences  which  have  been  dropped  off 
of  the  tree  are  insignificant.  This  is  called  hypothesis  merging.  (A  more 
thorough  treatment  of  the  foregoing  concepts  may  be  found  in  [6].) 

In  order  to  apply  the  above  methods,  we  must  have  some  means  of 
measuring  the  strength  of  the  various  hypotheses.  In  other  words,  we  must 
define  a  probability  measure  on  the  set  of  branches  of  the  tree.  The  par¬ 
ticular  definition  and  evaluation  of  a  probability  measure  depends  upon  the 
nature  of  the  tracking  system  and  the  types  of  measurements.  One  observa¬ 
tion  can  be  made  here,  however.  For  data  association  trees,  a  natural  de¬ 
finition  of  probability  is  one  that  is  monotonically  related  to  the  "close¬ 
ness"  of  a  measurement  to  the  predicted  value  of  a  data  track.  The  closer 
a  measurement  is  to  a  particular  prediction,  the  higher  the  probability 
that  the  measurement  is  associated  with  the  corresponding  data  track.  This 
use  of  predictors,  and  the  fact  that  we  are  using  sequential  processing. 
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suggests  the  use  of  Kalman  filters.  A  detailed  application  of  this  idea 
appears  in  Chapter  4. 

The  definition  of  this  "closeness  probability"  leads  to  the  conoept  of 
clustering.  (See  Fig  2-4).  Suppose  a  threshold  is  defined  such  that,  dur¬ 
ing  tree  expansion,  any  branch  whose  probability  lies  below  this  threshold 
is  pruned.  We  can  think  of  drawing  a  "gate"  around  the  predicted  value  of 
each  data  track.  The  probability  of  a  hypothesis  below  the  threshold  is 
then  equivalent  to  the  present  measurement  falling  outside  the  gate.  By 
considering  only  measurements  inside  the  gate,  we  have  effectively  pruned 
the  tree  of  hypotheses  with  probabilities  below  a  certain  level.  Suppose 
now  that  the  gates  of  various  data  tracks  do  not  overlap.  The  scan  can 
then  be  partitioned  into  a  number  of  subsets,  each  subset  containing  the 
measurements  falling  inside  the  gate  of  a  different  data  track.  Because 
the  possible  data  associations  of  each  of  these  subsets  are  mutually  ex¬ 
clusive,  we  can  break  up  the  data  association  tree  into  a  number  of  small¬ 
er,  independent  trees.  This  is  an  enormous  simplification,  for  the  sum  of 
all  of  the  hypotheses  in  the  smaller  trees  will  be  much  less  than  the 
number  of  hypotheses  in  the  tree  which  spawned  them.  For  fixed  resources, 
this  amounts  to  an  increase  in  the  number  of  hypotheses  that  can  be  con¬ 
sidered  simultaneously.  If  some  of  the  gates  do  overlap,  we  can  form  one 
tree  for  the  cluster  of  the  corresponding  data  tracks.  There  is  still  a 
simplification  in  this  case,  although  not  as  great.  The  concept  of  clus¬ 
tering  is  important  for  our  later  work. 


2J4  -  TRACK  HYPOTHESIS  REDUCTION  -  THE  DELAYED  il-SCAH  ALGORITHM 

The  data  association  trees  of  different  tracking  systems  are  usually 
quite  similar  in  structure.  The  construction  of  track  association  trees 
will  vary  widely  with  the  tracking  system.  In  some  cases,  the  track  asso¬ 
ciation  trees  are  degenerate,  and  target  state  estimates  are  a  trivial 
consequence  of  data  association.  For  instance,  take  the  case  of  a  single 
node,  single  target  tracker  that  uses  position  measurements.  Initial  tar¬ 
get  states  can  De  estimated  U3ing  only  two  measurements.  In  fact,  the  tar¬ 
get  state  can  be  used  iri  a  Kalman  filter  to  evaluate  probabilities  in  the 
tree.  This  formulation  is  precisely  the  "N-sean  algorithm”  developed  by 
Singer,  Sea,  and  Housewright  [93.  The  term  "N-scan"  refers  to  the  fact 
that  the  hypothesis  tree  is  depth-limited  to  U-scans,  and  hypothesis  merg¬ 
ing  is  done  over  the  past  Il-scans.  A  generalization  to  the  multi-target 
case  was  presented  by  Reid  [<ij.  Here  again,  target  states  are  used  for 
probability  calculations.  Reid  uses  clustering  in  his  algorithm  a3  well  as 
an  N-scan  approach. 

The  generalization  of  these  algorithms  to  a  multi-node  active  system, 
of  course,  requires  the  correlation  of  state  estimates  from  different 
nodes,  and  hence  requires  trees.  However,  assuming  that  node  estimates  arc 
independent  of  one  another,  and  that  noise  and  target  density  are  suffi¬ 
ciently  low,  ambiguities  will  resolve  themselvc-s  fairly  rapidly.  The  case 
is  different  for  passive  systems,  which  are  inherently  multi-node  in  na¬ 
ture.  Even  if  it  is  possible  to  estimate  target  state  from  the  measure¬ 
ments  of'  a  3ingle,  passive  node,  the  covariance  matrix  of  such  an  estimate 
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is  usually  so  largo  that  estimate  is  virtually  useless.  Correlation 
between  data  tracks  of  different  nodes  in  the  system  is  thus  necessary  to 
obtain  good  estimates. 

Correlation  of  data  tracks  requires  time,  and  so  reduction  of  track 
association  trees  is  a  rather  drawn  out  process.  It  would  be  desirable  to 
do  as  few  of  these  operations  ns  possible.  This  is  the  motivation  for  the 
following  scheme,  called  the  "delayed  N-scan  approach."  Instead  of  con¬ 
structing  track  association  trees  for  each  data  association  hypothesis 
under  consideration,  we  delay  the  construction  until  data  tracks  are  deter¬ 
mined  in  an  N-scan  algorithm  on  the  oata  association  tree.  In  other  words, 
a  measurement  is  not  used  in  updating  probabilities  of  track  hypotheses  un¬ 
til  the  scan  to  which  it  belongs  becomes  identified  in  the  data  association 
tree.  This  technique  effectively  separates  the  data  association  and  track 
association  processes.  The  resulting  benefits  are  the  same  as  in  employing 
clustering,  since  the  track  association  trees  were  originally  nested  in  the 
data  association  tree.  The  disadvantage,  of  course,  ip  that  the  incorpora¬ 
tion  of  measurements  into  state  estimates  is  delayed  by  the  N-scan  algo¬ 
rithm  on  the  data  association  tree.  It  is  the  price  paid  for  simplifica¬ 
tion  of  the  problem.  This  procedure  will  be  used  in  our  application  in 
Chapter  jj. 

2.5  -  OTHER  WORK  IN  MULTITARGET  TRACKING 

Before  leaving  this  chapter,  we  should  summarize  previous  work  in  inul- 
titar&et  tracking  and  compare  it  to  the  approach  taken  here.  An  excellent 
survey  of  multitarget  tracking  methods  can  be  found  in  a  paper  by  Bar- 
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All  of  the  work  in  the  literature  focuses  primarily  on  developing 
techniques  to  resolve  the  uncertainty  in  the  origin  of  received  measure¬ 
ments;  i.e.,  the  data  association  problem.  The  applications  of  these  tech¬ 
niques  generally  assume  that  target  state  estimates  can  be  obtained  direct¬ 
ly  from  given  sequences  of  measurements. 

The  approaches  to  the  problem  can  be  classified  as  Bayesian  and  non- 
Bayesian.  Non-Bayesian  approaches  do  not  take  into  account  a  priori  infor¬ 
mation.  The  early  work  of  Sittler  [10]  typifies  this  approach.  In  his  al¬ 
gorithm,  data  association  trees  are  formed  in  a  similar  manner  as  in  our 
algorithm,  although  this  is  not  explicitly  stated.  Kalman  filters  estimat¬ 
ing  target  states  are  initiated  and  updated  directly  by  the  given  measure¬ 
ments.  The  innovations  of  the  filter  of  each  branch  are  used  to  sequen¬ 
tially  compute  a  likelihood  function.  Target  tracks  whose  likelihoods  are 
below  a  certain  threshhold  value  are  then  rejected.  A  somewhat  different 
technique  developed  by  Morefield  L7J  minimizes  the  likelihood  function  by 
transforming  the  problem  into  an  integer  programming  problem.  This  pro¬ 
duces  the  most  likely  set  of  target  tracks  given  all  of  the  data,  but  is  a 
batch  processing  technique. 

Two  observations  can  be  made.  The  first  is  that  the  state  estimates 
and  covariances  are  conditioned  on  the  corresponding  data  hypotheses  being 
true.  However,  no  probability  value  is  obtained  for  the  data  hypotheses 
themselves.  This  Is  essentially  due  to  the  philosophy  of  the  non-Bayesian 
approach.  The  second  observation  is  that,  without  a  priori  information,  it 
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is  difficult  to  apply  these  algorithms  in  a  distributed  processing  system, 
because  the  data  associatibn  and  track  association  are  intertwined. 

The  Bayesian  approaches  to  th,e  multitarget  tracking  problem  attach  to 
each  measurement  a  probability  of  being  correct,  based  on  a  priori  informa¬ 
tion.  The  resulting  target  state  estimates  and  covariances  reflect  the  un¬ 
certainty  in  the  origin  of  the  measurements.  In  [1],  Bar-Shaiom  and  Tse 
deal  with  the  single  target  case.  Assuming  a  target  state  has  already  been 
initiated,  a  gate  can  be  formed  around  the  estimated  target  state.  For 
each  measurement  that  falls  within  the  gate,  a  probability  of  being  correct 
is  computed  based  on  how  close  the  measurement  is  to  the  estimate.  These 
measurements,  weighted  by  their  probabilities,  are  then  used  to  update  the 
estimate.  The  resulting  filter  is  called  the  probabilistic  data  associa¬ 
tion  filter  (PDAF).  This  approach  is  target  oriented  in  nature.  A3  such, 
it  is  difficult  to  incorporate  initialization  of  target  tracks  into  the 
scheme.  The  PDAF  is  extended  to  the  multiturget  case  in  [23.  However,  the 
scheme  for  computing  probabilities  is  very  complicated.  In  addition,  the 
application  of  these  ideas  to  a  distributed  processing  system  would  in¬ 
crease  enormously  the  amount  of  computation  required.  Because  of  this,  and 
our  desire  to  include  track  initiation,  we  have  rejected  this  approach. 

The  work  of  Dinger  et.  al.  i-J J  and  Reid  [53  has  already  been  men¬ 
tioned.  We  have  described  the  delayed  N-scan  algorithm  as  a  generalization 
of  Reid.  It  is  probably  more  accurate  to  characterize  it  as  a  combination 
of  the  Bayesian  and  non-Baye3ian  approaches.  In  the  sequel,  we  correlate 
data  tracks  in  the  track  association  trees  by  computing  a  liklihood  func¬ 
tion  for  each  track  hypothesis.  The  target  state  estimation  is  conditioned 


on  a  particular  data  hypothesis  (or  hypotheses)  being  true.  This  is  exact¬ 
ly  the  non-Bayesian  approach.  However,  the  data  hypotheses  upon  which 
these  target  tracks  are  conditioned  are  obtained  by  a  "pre-processor"  which 
employs  a  Bayesian  N-scan  algorithm  similar  to  Reid's.  Thus,  although  the 
target  state  estimates  do  not  reflect  measurement  origin  uncertainty,  they 
are  based  on  hypotheses  that  have  been  singled  out  in  a  Bayesian  weeding 
process . 

2.6  -  SUMMARY 

In  this  chapter ,  we  have  defined  some  basic  terms  and  constructed  a 
model  lor  the  solution  of  the  trucking  problem.  Ue  have  considered,  in  a 
general  setting,  some  useful  techniques  for  rendering  the  model  amenable  to 
practical  implementation,  and  we  have  seen  how  some  previous  work  in  track¬ 
ing  fits  into  this  model.  In  the  next  chanter,  we  present  an  application 
of  these  ideas  to  a  specific  trucking  system. 


CHAPTER  3 


DISTRIBUTED  SENSOR  NETWORKS 

The  procedures  developed  in  the  last  chapter  ere  quite  general  in  na¬ 
ture.  In  order  to  demonstrate  their  utility,  we  need  to  have  a  specific 
application.  In  this  chapter,  we  will  describe  a  general  tracking  system 
called  a  Distributed  Sensor  Network.  This  tracking  system  will  provide  the 
motivation  for  the  example  to  which  we  shall  apply  our  algorithm. 

In  principle,  the  more  information  a  tracking  system  can  obtain  from 
the  environment,  the  better  it  will  be  able  to  track  targets.  Active  sys¬ 
tems  are  almost  always  used  in  situations  where  maximum  target  information 
is  the  only  criteria  or  has  absolute  priority.  Host  active  systems,  such 
as  radar,  send  a  signal  into  the  environment  and  receive  reflected  versions 
of  those  signals,  along  with  clutter  from  the  environment.  This  operation 
of  reflection  allows  an  estimation  of  the  time  delay  between  transmission 
and  reception.  Since  the  speed  of  propagation  of  the  signal  is  assumed 
known,  this  is  equivalent  to  a  range  estimate.  Passive  systems,  on  the 
other  hand,  do  not  have  this  information,  and  hence  must  contain  more  sen¬ 
sors  and  nodes  than  active  systems  to  obtain  equivalent  information. 

In  many  cases,  however,  maximum  target  information  is  not  the  sole 
criteria,  nor  does  it  always  have  the  highest  priority.  Active  systems 
have  the  property  that  potential  targets  in  the  environment  may  be  able  to 
detect  the  radiated  signals.  In  many  cases,  especially  in  military  appli¬ 
cations,  this  detectability  is  a  major  drawback,  and  the  design  of  passive 
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systems  is  desirable. 

DSN  (Distributed  Sensor  Network)  is  a  research  project  of  MiT  Lincoln 
Laboratory.  DSN  is  a  inuiti-node  surveillance  system  designed  to  detect, 
locate,  track,  and  identify  low-flying  aircraft.  Although  in  actual  prac¬ 
tice  the  system  may  be  a  mixed  one,  with  both  active  and  passive  sensors, 
it  turns  out  to  be  more  fruitful  to  study  a  strictly  passive  system,  since 
the  distributed  processing  and  control  problems  of  large  act.  *e  systems  oc¬ 
cur  in  much  smaller  passive  systems.  In  addition,  in  reference  to  detecta¬ 
bility,  it  would  be  desirable  for  the  system  to  be  capable  of  carrying  out 
its  functions  using  only  passive  sensors.  Hence,  we  will  view  the  DSN  as 
having  a  solely  passive  capability. 

The  new  information  input  to  the  DSN  are  acoustic  signals.  Each  node 
of  the  DSN  contains  an  array  of  acoustic  sensors  (probably  high  quality  mi¬ 
crophones).  Each  sensor  samples  the  incoming  signals  at  a  specified  rate 
to  obtain  a  digital  sample  set.  After  predetermined  intervals  of  time,  the 
sampled  signals  from  all  sensors  are  input  to  the  signal  processing  com¬ 
ponent  of  the  node.  Essentially,  this  component  uses  high  resolution  fre¬ 
quency  wavenumber  analysis  to  detect  phase  differences  in  the  signals, 
thereby  obtaining  a  direction.  Tiie  output  can  be  viewed  as  a  curve  plot¬ 
ting  received  signal  power  vs  azimuth.  An  example  plot  is  shown  in  Figure 

3-1. 


These  power-azimuth  curves,  produced  at  given  instants  of  time,  are 
the  input  data  to  the  trucking  system.  To  keep  things  simple  for  our 
models  to  come,  we  will  use  only  the  azimuthal  information .  The  power  in- 
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formation  could  be  useful  iri  resolving  ambiguities;  however,  we  shall  nob 
consider  this  role. 

We  are  thus  assuming  that  each  node  of  the  D31I  has  available  a  set  of 
azimuths  and  variances  of  those  azimuths  at  disorete  instants  of  time.  The 
tracking  problem  is  to  combine  the  azimuths  from  the  nodes  to  produce  tar¬ 
get  tracks.  We  should  keep  in  mind  that  the  nodes  arc  geographically 
dispersed  and  that  the  velocity  of  signal  propagation  is  comparable  to  the 
speed  of  potential  targets.  In  addition,  the  measurement  times  at  dif¬ 
ferent  nodes  are  not  necessarily  synchronous.  However,  for  simplicity,  we 
shall  ignore  this  last  observation. 
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In  a  large  scale  DSN,  it  is  obvious  that  a  distributed  processing 
scheme  is  much  more  desirable  than  a  central  processing  one.  First  all, 
the  large  amounts  of  information  dictate  that  it  be  handled  in  pieces  in 
order  to  sort  it  out.  Second,  communication  is  generally  much  slower  than 
computation,  arid  so  for  time  efficiency,  the  computational  load  should  be 
distributed  as  much  as  possible.  Communications  also  radiate  power,  which 
is  undesirable  if  detectability  is  an  issue. 


Thus,  it  seerns  logical  to  do  as  much  of  the  processing  at  individual 
nodes  as  possible.  There  is  not  enough  information  gathered  by  the  node  to 
produce  independently  a  reliable  target  track.  However,  it  i3  possible  for 
the  node  to  perform  its  data  association  independently  of  other  nodes. 
This  is  why  we  separated  out  data  association  and  track  association  in  our 
algorithm.  It  fits  very  naturally  into  the  problem  of  distributing  compu¬ 
tational  load  in  a  tracking  system. 
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To  produce  the  actual  target  tracks,  wo  can  correlate  date  tracks  of 
pairs  of  nodes  to  produce  target  tracks.  Target  trucks  cun  then  be  refined 
at  a  higher  level  by  comparing  various  two-node  results.  This  scheme  sets 
up  a  hierarchy  similar  to  multi-site  radar,  with  the  exception  that  the  set 
of  two-node  target  tracks  are  not  always  independent. 

The  two-node  tracking  problem  is  what  we  intend  to  study.  We  will  as¬ 
sume  that  our  tracking  system  consists  of  two  nodes  each  independently  re¬ 
ceiving  azimuth  measurements.  For  simplicity,  we  will  assume  that  the  en¬ 
vironment  is  two  dimensional  and  that  nodes  and  targets  are  dimensionless 
points.  All  noise  in  the  measurements  is  assumed  to  be  Guassian  and  white. 

With  the  specification  of  our  example,  we  are  now  ready  to  explore 
mathematical  details.  The  first  order  of  business  is  to  discuss  a  defini¬ 
tion  of  probability  for  our  hypothesis  trees.  This  is  the  subject  of  the 


next  chapter. 
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CHAPfEii  4 

THE  CALCULATION  OF  HYPOTHESIS  PROBAUILITIES 

In  order  to  apply  our  algorithm  to  the  two-node  system,  vie  first  need 
to  define  a  probability  measure  on  both  the  data  association  trees  and  the 
track  association  trees.  With  this  we  can  evaluate  the  strength  of  various 
hypotheses  and  eliminate  unlikely  ones.  In  this  chapter,  ue  develop  a 
theoretical  basis  for  calculating  the  probabilities. 

4.2  -  DATA  ASSOCIATION  TREES 

We  first  consider  the  data  association  tree  at  a  single  node.  Sup¬ 
pose,  for  the  moment,  that  an  estimate  of  the  target  state  is  available  at 
a  measurement  time  t.  We  can  then  form  a  prediction  of  the  incoming  meas¬ 
urement  based  on  this  target  state  estimate.  Intuitively,  we  would  expect 
that  the  closer  the  measurement  is  to  the  predicted  value,  the  more  likely 
the  measurement  is  associated  with  the  target.  As  is  well  known,  we  can 
sequentially  form  target  3tate  estimates  as  data  arrives  through  the  use  of 
Kalman  filters. 

In  a  more  general  setting,  we  assume  that  the  measurements  are  the 
output  of  a  linear  system  plus  un  additive  noise  term.  As  is  well  known, 
we  can  formulate  a  state  variable  description  of  the  system  in  many  dif¬ 
ferent  ways  by  choosing  different  definitions  of  the  state  variables  of  the 
system.  For  the  moment,  let  us  assume  we  have  settled  on  a  particular  de¬ 
finition  for  the  state  of  the  system.  We  may  write  the  (linear)  state 
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equation  and  measurement  equation  us 

x(k+1)  =  F(k)x(k)  +  G(k)w(k) 
z(k)  s  H(k)jc(k)  +  _v< k)  <4 . 1  > 

x,Ck)  is  the  state  of  the  system  at  time  k  and  _z(k)  is  the  measurement  at 
time  k.  We  assume  that  w  and  v  are  independent  Gaussian  white  noise  se¬ 
quences  with 

E[w(k)wT(k)]  =  Q(k) 

E[v(k)vT(k)J  s  R(k)  (4.2) 

Define  jUkjl)  to  be  the  linear  least  squares  estimates  of  jc(k)  given  data 
up  to  time  1  and  P ( k i 1 )  to  be  the  covariance  matrix  of  thi3  estimate.  We 
can  obtain  jc(kik)  and  P(kik)  through  the  use  of  the  discrete  time  Kalman 
filter  equations  as  follows: 

jf(kik)  =  jc(k|  k-1 )  +  K(k)[z(k)  -  H(k)X(ki k-1 ) j 
£(k+1ik)  s  F(k)8(k|k) 

K(k)  =  P(k!  k-1  )HT(k)  L1I (k) P (k !  k— 1  )HT(k)  +  K(k)]"1 
P(kik)  =  [I  -  K(k)H(k)]P(k|k-1 )[I  -  K(k)IKk)]  +  K(k)R(k)KT(k) 
P(k+1 i k)  =  F(k)P(kik)FT(k)  +  G(k)Q(k)GT(k) 

(4.3) 

Using  these  equations,  we  can  sequentially  update  the  state  estimate 
as  measurements  arrive.  One  quantity  that  will  be  of  interest  to  us  is  the 
so-called  innovations  sequence 


</k  =  z(k)  -  H(k)jUkik-l) 


,ij wo*  •; rc.i yki  n  '  , — "»n,  ittii  u,>  in- 


It  can  be  shown  (e.g.  L 1 1 3 )  that  is  a  Guussian  white  noise  sequence 
given  our  previous  assumptions  with 

EtC^J]  =  Bk  =  H(k)P(kik-1  )HT(k)  +  R(k)  (4.4) 

V/e  are  now  ready  to  derive  our  result  concerning  the  calculation  of 
probabilities  of  hypotheses.  Uur  deviation  follows  closely  the  work  of 
Reid  Ld].  Although  Reid  assumes  that  jc(k)  is  the  actual  target  state,  the 
results  are  valid  for  our  more  general  definition  of  x.(k). 


Let 


2(k)  =  {z  (k)  ,  in=1 ,2 . M,  } 

•’ll!  K 

denote  the  measurements  received  at  time  k  and 

Zk  =  LZ(  1 )  ,Z(2) . Z(k) } 


denote  the  cumulative  set  of  measurements  through  time  k.  Also,  define 

V  k  =  tr  k  .  i=1 ,2 . ik> 

*  '  k 

to  be  the  cumulative  set  of  hypotheses  just  after  time  k.  Each  r  ^ 

corresponds  to  s  brunch  on  the  hypothesis  tree. 


Now,  define 


Pk  =  P(r  k  i  Zk) 


k  k 

that  is,  is  the  probability  at  time  k  of  the  branch  T  ^  of  the  hypothesis 

tree  given  the  data  through  time  k  (Zl<).  In  actuality,  this  is  equivalent 

k-1 

to  the  conditioned  Joint  probability  of  the  prior  hypothesis  V  f  and  the 
data  hypothesis  for  the  current  measurement  "el  i}i  .  Dropping  the  depen¬ 
dence  on  pa3t  data  for  notutional  simplicity,  we  ctin  use  Dayes1  equation  to 


write  the  relationship 


pi  s  p(rs’1,Sih  1  z(k)>  e  ip(z(k)  1  r£"1,^)  x  p%  1  rg'1) 


(4.5) 


The  term  rj  is  a  normalizing  factor  given  by 

k 
B 


{]  =  r  t  p<z(k)  i  r  k-1  ,ip. )  x  p<*.  !  r  J"1 )  x  p(r  jj”1 ) 

g  h  B  r»  n  g  B 


If  we  can  find  expressions  for  the  first  two  terms  on  the  right  hand  side 
of  (4.5),  we  will  have  a  recursive  relationship  for  calculating  probabili¬ 
ties  . 


The  first  term  is  the  probability  density  function  of  the  current  set 

k-1 

of  measurements  Z(k)  given  the  prior  hypothesis  V  r  ,  and  the  current  data 

B 

hypothesis  Assuming  that  each  measurement  z^(k)  in  Z(k)  is  condition¬ 

ally  independent,  we  have 

Mk 

p(zik)  i  r  k“\i|ih)  =  n.  p(z  (k)  i 

8  n  ns  1  b  11  (4.6) 

Suppose  that  41^  assigns  z^OO  as  either  a  false  alarm  or  a  new  data 
track.  In  either  case,  there  is  no  a  priori  information  to  determine 
whether  one  set  of  possible  measurements  is  more  likely  than  another. 
Hence,  in  these  two  cases,  we  will  assume  that  i'c  appears  as  a  uniform  dis¬ 
tribution 


P  ( z 

Vk> 


(4.7) 


V  is  the  volume  (or  area)  of  the  part  of  the  environment  covered  by  the 
node . 


If  neither  case  above  holds,  then  d.  assigns  z  (k)  to  cither  a  previ- 

Mi  m 


ously  established  aata  track  (confirmed  track)  or  to  a  data  track  whose  ex- 


k-1 

istence  is  implied  by  T  (tentative  track).  Tentative  traek3  are  ini- 

O  7 

tiated  when  a  new  data  track  is  hypothesized  in  the  course  of  expansion  of 
the  data  association  tree.  A  tentative  track  becomes  confirmed  if  it  still 
exists  when  the  scan  in  which  it  is  initiated  is  identified. 


If  a  Kalman  filter  is  running  on  the  assigned  data  track,  we  have 
available  a  current  estimate  8(k|k)a8k.  Since  this  estimate  does  not 
depend  on  the  current  measurement,  we  have 

p<Vk)  1  rf1)  a  P^Ckj-HCk)^  i  rk”\inh}  <n.8) 

The  last  is  the  conditional  density  function  of  the  innovations  of  the 
filter,  which  has  a  normal  distribution.  Hence 

PU^k)  i  =  M[<^n(k)-H(k)Jk,Bk]  (4.9) 

with 


Bk  a  H(k)P(k!k)HT(k)  +  R(k) 


N  (.(/,?)  = 
where  n 


exp(-^W) 


is 


sion  of  the  innovation  vector  x. 


(4.10) 


The  second  term  on  the  right  hand  side  of  (4.5)  is  the  probability  of 
the  current  data  association  hypothesis  given  the  prior  hypothesis  T  k_1 


8 


has  three  items  of  information: 


a)  Number  -  the  number  of  measurements  associated  with  prior  data 
tracks  (N^,(h)),  false  alarms  (Np^,(h))  and  new  data  tracks  (H^(h)). 

b)  Configuration  -  the  partition  of  the  set  Z(k)  into  three  subsets 
corresponding  to  prior  targets,  false  targets,  and  new  targets. 

c)  Assignment  -  the  assignment  of  each  measurement  associated  with 
prior  data  tracks  to  the  specific  source. 
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In  addition,  the  prior  hypothesis  gives  N^(g),  the  total  number  of 
data  tracks,  confirmed  and  tentative,  implied  by  that  hypothesis.  Thus  if 
N^N,^ ,  yh  implicitly  assumes  that  measurements  were  not  received  from 
some  of  the  prior  data  tracks. 

k-1 

To  find  the  probability  of  the  numbers  N^,  hp^,  arid  given  I*  ,  we 
make  the  following  assumptions.  First,  the  probability  that  a  target  will 
generate  a  measurement  which  is  actually  received  by  the  node  is  a  oonstant 
Pp  (called  the  probability  of  detection).  In  other  words,  the  reception  of 
a  measurement  can  be  described  in  probabilistic  terms  as  a  Bernoulli  trial. 
If  we  further  assume  that  the  detectability  of  each  target  is  independent 
of  the  others,  and  that  each  target  can  only  generate  one  measurement  in  a 
given  scan,  then  we  see  that  Np_.  is  a  Bernoulli  prooe33  and  its  distribu¬ 
tion  is  binomial. 

Second,  we  assume  that  the  number  of  false  alarms  follows  a  Poisson 
distribution.  This  is  ari  assumption  often  made  in  tracking  problems.  It 
makes  intuitive  sense,  because  while  the  appearance  of  false  alarms  is  ran¬ 
dom,  in  many  cases  they  have  a  constant  average  rote  of  appearance  over 
reasonably  lengthy  periods  of  time.  We  will  also  assume  the  number  of  new 
targets  follows  a  Poisson  distribution.  The  assumption  is  much  harder  to 
justify.  In  order  to  make  it  reasonable  in  actual  practice,  the  average 
rate  of  new  target  appearances  must  be  adjusted  much  more  often  than  the 
rate  for  false  alarms. 

With  these  assumptions,  we  have 
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p(NDT,N:jT,r,FT  1  r  a  )  \nm)  (1“V 


X  FHFT(^FTV)FN1,JT(fiNTV) 


(4.11) 


where 


Pp  =  probability  of  detection 
pFT  s  density  of  false  alarms 
PWT  =  density  of  previously  unknown 


targets  that  have  been  detected 


vn  -X 

F  (X)  = 
n  n  1 


Now,  the  total  number  of  measurements  is 


Hk  =  MDT  +  ;JFT  +  UHT 


The  number  of  different  partitions  of  the  M  measurements,  given  the 


numbers  ,  Np^,  and  N^T  is 


M  V1JDT 


ndt/V  nft 


Assuming  that  each  configuration  is  equally  likely,  the  conditional  proba¬ 


bility  of  a  specific  configuration  is 


P(Configuration  !  f,DT •  NpT , NJJT >  = 


(4.12) 


'ot/JW 


Given  the  configuration,  the  number  of  possible  assignments  of  the  N,. 


AJT61  ■ 


measurements  to  the  IJ  prior  data  tracks  is 

lul 


Assuming  that  each  suoh  assignment  is  equally  likely,  we  have 

(N  -M-..,)! 

P(Assignment  i  Configuration)  =  — — 

htgt  } 


(4.13) 


The  joint  conditional  probability  of  NDT ,  hpT ,  ,  the  configuration,  and 
the  assignment  is  the  product  of  (4.11),  (4.12),  (4.13),  and  is  also  the 
conditional  probability  of  i|i.  .  Thus  we  get 


U— 1  ^MT  ^  Nrvr  (  NTf,T-N nT  ^ 

p(ei  i  rK  ')  s  — — ill-  x  p  DTd~p  )  TC,T  DT 

KVh  '  1  g  1  Mk!  X  *£>  u  *V 


FHFT(fjFTV)FtJNT(^NTV) 


(4.14) 


If  we  now  substitute  (4.6),  (4.7),  (4.9)  and  (4.14)  into  (4.5)  we  obtain 


k  _  1  HFT!tW  NDT  .  (fJTGT“l’iDT) 

Fi  “  t|  Mk!  X  u  (,-V 


x  FNft(^FTV)F«n't(^NTV) 


x  IJ  N[z  (k)-H(k)3k  ,  Bk]  -1-  ■  ■ 
m.1  ",n  “  vUFT+NNT 


(4.15) 


The  measurements  have  been  implicitly  ordered  in  the  above  so  that  the 
first  NpT  of  them  correspond  to  those  assigned  to  prior  data  tracks  by 
Substituting  into  (4.15),  simplifying,  and  incorporating  constants 
in  rj,  we  finally  get 


p*  .  lpNDT(1  p  (NTGT"NDT)NFT"NT 
Pi  -  r)PD  °-PD)  ^FT  Mut 

ndt 

X  n  »tz  (k)-il (k) 5?k  ,  bkJ  P^1 

rn=  1  -m  •  m  4 


(4.16) 


Hote  that  this  is  independent  of  V. 


As  mentioned  by  Raid,  this  equation  is  easily  implemented  in  the 
framework  of  tree  expansion.  Prior  to  expansion,  all  prior  probabilities 

ntgt 

are  multiplied  by  (1-Pp)  .  Then,  as  each  branch  is  created  for  a 

specific  measurement  assignment,  we  multiply  the  probability  of  the  prior 
hypothesis  by  either  pfT,  p^,  or  P^ltz^CkJ-HCk) ,  B^j/U-Pp),  depending 
on  the  assignment.  The  probabilities  can  be  normalized  after  tree  expan¬ 
sion,  although  this  is  not  strictly  necessary  since  only  the  relative  pro¬ 
babilities  are  important. 


-ihPft 
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for  false  alarms 
for  new  data  tracks 


§ln2ir  ♦  ^ln|B|  -  lnPD  +  lnd-Pp) 


for  prior  data  tracks 


We  see,  then,  that  we  need  to  specify  three  items  of  information  in 
order  to  use  this  definition  on  the  data  association  trees: 

a)  The  state  variable  representation  (4,1) 

b)  The  false  alarm  density  p^ 

c)  The  new  target  density  p^ 


^’  \vr  PfP-  "SU'W r-w  -  .T'TvmrT1 
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We  next  consider  track  association  trees. 

4.2  -  TRACK  ASSOCIATION  TREES 

Track  association  trees  behave  somewhat  differently  than  data  associa¬ 
tion  trees.  Here  we  hypothesize  various  combinations  of  data  tracks  from 
different  nodes  to  produce  possible  target  tracks.  Assuming  that  each  such 
combination  produces  a  unique  target  track,  we  see  that  these  trees  hove 
constant  depth.  When  a  new  data  track  appears,  new  nodes  ore  created  by 
correlating  the  new  track  with  existing  tracks  from  other  nodes.  The  pro¬ 
cess  is  a  breadth  expansion,  rather  than  depth  because  the  new  hypotheses 
are  not  conditioned  on  the  hypotheses  already  in  existence.  There  is  no 
need  to  consider  false  alarms  or  new  targets  at  this  level;  these  have  al¬ 
ready  been  determined  at  the  data  association  level.  Wo  thus  need  only  a 
method  of  distinguishing  those  combinations  of  data  tracks  that  correspond 
to  real  targets.  Those  that  do  not  correspond  to  real  targets  are  termed 
ghosts . 

The  problem  in  setting  up  a  general  definition  of  probability  here  is 
that  the  probability  of  each  node  must  be  evaluated  over  a  time  interval. 
Tnis  differs  considerably  from  the  nodes  in  data  association  trees,  whose 
probabilities  aru  evaluated  at  points  in  time.  If  new  nodes  are  added  to 
an  existing  track  association  tree,  then  the  probabilities  of  the  new  nodes 
will  be  evaluated  over  different  time  intervals  than  the  older  ones.  The 
question  then  becomes  one  of  comparing  these  probabilities. 


Instead  of  computing  probabilities,  we  will  compute  a  likelihood  func¬ 
tion  for  each  track  hypothesis.  Suppose  we  are  using  the  data  tracks  from 
n  nodes  to  form  a  track  hypothesis.  Define  ^(k)  to  be  the  augmented  meas¬ 
urement  vector,  consisting  of  the  measurement  vectors  of  the  n  nodes.  We 
can  model  the  problem  as  follows: 

x(k+1)  =  A(k)x(k>  +  B(k)w(k) 
z(k,0)  =  C(k)jc(k)  +  ^(k) 

k)  here  is  the  target  state.  The  parameter  6  denotes  the  dependence  of  z 
on  the  track  hypotheses.  To  find  the  most  likely  hypothesis,  wo  need  to 
maximize  the  likelihood  function  for  0. 

The  problem  in  this  form  resembles  the  multiple  model  identification 
problem  (see  Van  Trees  till).  Using  the  model,  we  can  obtain  a  recursive 
relationship  for  the  log  likelihood  function  of  6.  The  result  is 
L(0,k+1 )  =  L(0,k)  +  (k,0)£T(k,0) 

where  C^(k,0)  is  the  innovations  at  time  k  obtained  from  a  Kalman  filter  as¬ 
sociated  with  the  track  hypothesis  identified  with  0,  and  B(k,0)  is  its  co- 
variance.  With  this  relationship,  we  have  a  means  of  choosing  the  most 
likely  track  hypothesis. 

The  initialization  of  a  track  hypothesis  occur:-,  when  a  new  data  track 
appears  out  of  the  data  association  process  at  some  node.  The  initializa¬ 
tion  becomes  somewhat  complex  when  the  measurements  are  time  delayed.  In 
some  cases,  it  may  not  be  possible  to  initialize  the  correct  track  hy¬ 
pothesis.  For  instance,  consider  the  situation  in  which  a  target  is  closer 
t©  a  -n  >de  A  than  .'mother  node  B.  Because  of  the  signal  propagation  ’ay, 
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the  measurements  ut  A  will  be  more  recent  then  those  at  B.  Presumably, 
node  A  will  detect  the  target  first.  If  we  attempt  to  correlate  the  new 
data  track  with  existing  data  tracks  ut  D,  we  see  that  the  correct  correla¬ 
tion  will  not  be  formed,  since  13  has  not  yet  detected  the  target. 

With  this  in  mind,  we  present  the  following  scheme  for  the  construc¬ 
tion  of  track  association  trees.  When  a  node  initiates  a  new  duta  track 
tt ..  ough  its  data  association  process,  a  new  track  association  tree  is  also 
created.  We  attempt  to  initialise  correlations  with  existing  data  tracks 
from  the  other  nodes  for  a  given  length  of  time.  After  this  period,  no 
more  attempts  at  initialization  are  made.  In  this  manner,  the  correct 
track  hypothesis  should  be  formulated  with  the  last  node  to  detect  the  tar¬ 
get.  In  order  to  compare  hypotheses  on  different  truck  association  trees, 
we  shall  use  time  averaged  values  of  likelihoods  computed  in  the  manner 
described  above;  i.e.  we  shall  use  L(k,d')/At,  where  At  is  the 

The  utility  of  the  above  method  will,  of  course,  depend  on  the  track¬ 
ing  system.  In  some  cases,  it  mny  turn  out  that  differentiating  between 
ghosts  and  targets  is  virtually  impossible.  We  then  must  either  bring  to 
bear  other  information  in  the  system  to  distinguish  the  real  targets,  or  we 
must  continue  to  track  the  ghosts  as  targets.  These  issues  will  appear  in 
our  discussion  in  the  next  chapter. 

J4.2  -  SUMMARY 


In  this  chapter,  we  have  presented  a  formulation  for  the  calculation 
of  probabilities  of  hypotheses  ori  the  data  association  trees  and  the  track 


association  trees.  Armed  with  this  knowledge,  we  can  now  proceed  with  the 
application  of  the  algorithm  to  the  two-node  system.  The  necessary 
mathematics  are  developed  in  the  next  chapter. 
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CHAPVER  5 

STATE  SPACE  MODELS  FOR  THE  TWO  NODE  SYSTEM 

In  order  to  apply  the  formulation  presented  in  Chapter  we  need  to 
develop  state  variable  representations  to  implement  the  neoessary  Kalman 
filters.  Our  goal  in  this  chapter  is  to  develop  the  appropriate  state 
space  models  for  both  the  data  association  process  and  the  track  associa¬ 
tion  prooess.  Although  we  wish  to  obtain  filters  that  perform  reasonably 
well,  we  do  not  undertake  a  thorough  examination  of  these  state  space 
models.  As  a  consequence,  some  important  issues  are  left  unexplored.  How¬ 
ever,  our  main  purpose  is  to  demonstrate  the  tracking  algorithm  developed 
earlier;  we  do  not  propose  to  derive  the  optimum  tracking  filters  for  the 
specific  tracking  system  under  consideration.  Hence,  we  content  ourselves 
with  n  less  detailed,  out  adequate,  analysis  focused  on  the  data  associa¬ 
tion  aspect  of  the  problem. 

5  .1  -  STATE  SFACE  MODEL  :  THl_  DATA  ASSOCIATION  FILTER 

The  most  natural  way  to  define  a  state  variable  for  this  problem  is 
the  target  state  as  defined  in  chapter  2.  It  is  easy  to  see  how  the  meas¬ 
urements  are  related  to  the  target  state.  Iri  Figure  5-1,  16  is  the  noise¬ 
less  measurement  (the  acoustic  azimuth)  at  time  k.  because  the  speed  of 
sound  is  finite,  this  means  that  6  corresponds  to  a  target  position  at  time 
t  sometime  in  the  past.  (This  is  called  the  acoustic  target  position  to 
distinguish  it  from  the  true  target  position  at  time  k) .  The  relationship 
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between  k  end  t  is 


k  a  t  ♦  £  (5.1) 

where  8  is  the  range  to  the  acoustic  target  position  (  the  acoustic  range) 
and  c  is  the  speed  of  sound,  liow,  let  |3  be  the  acoustic  azimuth  (at  time 
kg)  when  the  acoustic  target  position  is  at  the  closest  point  of  approaoh 
(CPA).  The  time  the  target  is  aotually  at  CPA  is  tQ.  Lot  y  be  the  angle 
around  CPA  -  that  is  the  angle  between  the  line  to  the  apparent  target  po¬ 
sition  P  and  the  line  to  the  CPA.  We  see  that 

-~  <  y  <  |  <5.2) 

Let  us  assume  that  tne  velocity  v  and  the  heading  is  constant.  The  dis¬ 
tance  between  P  and  CPA  is 

d  =  v(t-tQ)  (5.3) 

A  negative  distant-  implies  the  target  is  approaching  CPA.  We  also  have 

d  =  rtany  (5.4) 

where  r  is  the  distance  to  CPA.  Thus 

v ( t— tQ )  =  rtany  (5.5) 

or,  substituting  for  t  and  tQ, 


v(k  -  - 
c 

tany  + 


+  — )  =  rtany 

-  i 


But  R  =  rsecy,  and  so 

tany  +  -secy  =  — (k-kn)  ♦  — 
c  rue 

Solving  for  y,  and  substituting  into  (5.2),  we  finally  obtain 
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t  =  p  ♦ 

The  quantities  kn 

* 

Hence ,  we  could 
model . 


>-*0>  1  ♦  ™ +  Vk-“o>2  +  fr(k-ko) 

_  r 

1--^ 

2 
c 

■p,  andp  provide  full  knowledge  of  the  target  state, 
choose  these  to  be  the  state  variables  for  the  system 


There  are  two  major  difficulties  in  using  a  state  model  in  this  form. 

First  of  all,  it  is  difficult  to  obtain  an  estimate  of  p  until  past  CPA. 

Because  the  governing  equation  (5.6)  is  nonlinear,  we  must  resort  to 

linearized  version  of  the  Kalman  filter  (called  the  extended  Kalman  filter 

or  EKF).  The  quantity  p  is  crucial  for  uri  accurate  linearization.  The 
» 

best  way  to  obtain  an  estimate  of  p  is  to  set  up  a  bank  of  filters,  each 
conditioned  on  a  different  value  of  p.  Based  on  the  performance  of  these 
filters,  we  would  then  apply  some  sort  of  decision  criteria  to  selsect  the 
best  filter.  This  has  been  done  by  Hebbert  [5!i  for  the  case  of  a  single 
target,  no  system  noise,  and  an  infinite  signal  propagation  speed.  In  our 
case,  there  will  already  be  several  Kalman  filters  running,  and  replacing 
each  of  these  with  a  bank  of  filters  increases  tho  memory  and  computation 
enormously. 

The  second  major  difficulty  in  using  the  full  target  state  in  our 
model  is  the  resulting  weak  observability  of  the  system.  Theoreticd’l  we 
can  obtain  information  about  every  state  variable  from  any  given  measure¬ 
ment.  However,  we  get  a  relatively  large  amount  of  information  about  a 
particular  state  variable  only  at  the  expense  of  the  other  3tate  variables. 
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Essentially,  the  azimuth  6  gives  one  position  coordinate  in  a  two  dimen¬ 
sional  system.  We  would  expect,  then,  that  we  could  not  obtain  substantial 
reductions  in  the  uncertainty  of  both  position  coordinates  is  shown  clearly 
in  Figure  5-2 (a-d).  x  is  defined  as  the  position  coordinate  parallel  to  a 
given  target's  trajectory,  and  y  to  be  position  coordinate  perpendicular  to 
the  trajeotory.  v  and  v  are  the  velocities  in  the  x  and  y  directions 

a  y 

respectively.  Figures  5-2 (a-d)  display  the  square  root  of  the  variances  of 

x,  y,  v  ,  and  v  as  a  function  of  time  for  a  typical  target  trajectory, 
x  y 

The  variances  are  the  diagonal  elements  of  the  error  covariance  matrix  of 

an  EKF  tracking  x,  y,  v  ,  and  v  .  The  filter  was  always  linearized  about 

x  y 

the  exact  trajectory.  The  measurement  noise  and  the  system  noise  in  this 
filter  were  set  to  U  (i.e.,  Q=R=0) ,  and  the  init  ial  covariance  was  set  to  a 
very  large  value.  The  filter  thus  starts  oil’  with  essentially  no  a  priori 
information,  and  all  reductions  of  the  covariance  are  solely  due  to  the  in¬ 
coming  measurements.  The  covariance  matrix  obtained  is  known  as  the 
Cramer-Rao  lower  bound,  which  implies  that  no  tracking  scheme  can  attain 
mean  square  errors  less  than  those  shown  here*.  As  can  be  seen,  there  is  a 
rapid  reduction  in  the  variance  of  while  the  target  is  still  far  from 
CFA  (which  occurs  at  t=45  seconds) .  The  variance  of  vx  stays  relatively 
high  until  the  target  approaches  CPA.  Thus,  in  the  region  approaching  CPA, 
the  full  state  tracker  is  at  best  weakly  observable.  The  large  uncertainty 
in  al  least  one  coordinate  in  this  region  could  result  in  a  poor  lineariza¬ 
tion  and  filter  divergence.  Since  we  wish  to  perform  data  association  in 
this  area,  we  must  reject  the  full  target  state  as  a  model. 


We  thus  must  deline  our  state  variables  in  another  way.  One  method  is 
make  a  Taylor  series  approximation  for  6,  i.e., 

T2 

d(k+T)  =  d(k)  +  T«5'(k)  +  ^-d"(k)  ♦  ~i6'"(k)  +  ...  (5.7) 

Our  problem  is  to  determine  the  appropriate  number  of  terms  to  retain  in 
the  expansion.  A  typical  plot  of  the  acoustic  azimuth  vs.  time  is  shown  in 
Figure  5-3.  Obviously,  a  linear  approximation  is  not  very  good  over  the 
entire  range,  so  the  term  d"  should  et  least  be  retained.  On  the  other 
hand,  retention  of  too  many  terms  can  lead  to  a  rather  sluggish  filter  with 
a  relatively  high  average  error  covariance.  In  addition,  computation  and 
memory  requirements  increase  with  the  number  of  state  variables.  There¬ 
fore,  we  restrict  our  attention  to  two  possibilities:  the  three-state  vec- 


and  the  four-3tate  vector 


Our  system  model  is 


xA(k)  =  d'(k)  1 

U"(k) 


A’Ck) 


*U(k)  =  U"(k) 


x ( k+T )  =  Fx(k)  + 


(5. da) 


where 


50 


T|< 
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0 

0 

T 

1 

0 
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T 

1 
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T 

1 
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x(k)  =  xA(k) 


x(  k)  s  xBUO 


(5.8b) 

We  assume  that  w(k)  is  a  scalar,  zero-mean,  white  Guassian  noise  process 
with  variance  Q.  Noise  is  entered  only  through  the  highest  derivative  for 
simplicity.  The  measurement  equation  for  the  system  is 

z(k)  =  Hx(k)  +  v(k)  (5.9a) 


with 


_x(k)  =  £A(k) 
x(k)  =  £B(k) 


(5.9b) 


v(k)  is  a  scalar,  zero  -mean, white  Guassian  noise  pi$e^ess  independent  of 
w(k)  with  variance  R. 


The  above  model  does  have  a  serious  drawback,  however.  When  the  ratio 
is  large,  correspc  ding  to  small  range  and  high  velocity,  the  actual 
6-curve  approaches  a  step  function,  and  its  derivatives  approach  singulari¬ 
ty  functions.  The  model  is  very  poor  in  these  situations.  In  fact,  the 
filter  will  not  react  as  fast  as  the  changes  in  these  derivatives,  and  it 
may  well  lose  the  data  track.  Our  solution  to  this  difficulty  is  to  bypass 
it,  pointing  out  that  it  is  easily  detected  from  an  azimuth  history,  and 
hence  special  mechanisms,  which  will  not  be  developed  here,  can  be  invoked 
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to  deal  with  it.  Thus,  we  will  e  isider  only  tercet  trajectories  with  rea- 

y 

sonably  small  values  of  — . 

In  order  to  properly  choose  a  good  model,  we  test  a  number  of  possible 
values  of  Q  in  both  the  3-state  ana  the  4-state  filter  over  a  range  of  typ¬ 
ical  trajectories.  The  trajecorics  were  produced  by  varying  the  velooity 
and  the  distance  to  CPA,  the  two  controlling  parameters.  Figure  5-3  shows 
the  various  combinations  that  were  used.  For  each  hypothesized  model,  a 
set  of  data  from  each  target  was  produced  by  adding  measurement  noise  of  a 
standard  deviation  of  j°  to  the  exact  acoustic  azimuths.  The  Kalman  filter 
corresponding  to  the  model  was  run  separately  on  each  data  set,  and  an 
average  squared  error  (ase)  between  the  true  acoustic  azimuth  and  the 
filter  estimate  was  computed  for  each  track.  An  overall  aso  was  computed 
for  each  filter. 

The  results  of  the  above  Monte  Carlo  simulation  are  shown  in  Figure 
5-4.  It  turned  out  that  filter  runs  from  two  of  the  trajectories  produced 
average  squared  errors  that  biased  those  computed  from  other  trajectories. 
The  first  .  ible  shows  the  results  if  these  runs  were  retained.  The  second 
run  gives  ol  computations  if  these  runs  are  eliminated  from  consideration. 
As  can  be  seen,  the  model  with  the  minimum  squared  error  is  different  in 
each  table.  The  erl'brs  are  close,  however,  indicating  that  there  are 
several  models  that  will  give  nearly  the  same  performance.  V/e  choose  the 
model  indicated  by  the  second  table,  because  all  of  our  example  target  try- 

y 

jectories  will  have  —  <  .1.  Thus,  we  choose  the  three  state  filter  with 
QsO.OOl  us  our  data  association  filter. 
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One  lust  item  to  consider  is  t  e  initialization  procedure  for  the  Kal¬ 
man  filter  based  on  (5. 8)  end  (5.y).  Assuming  the  model  is  reasonably  ob¬ 
servable,  the  filter  will  be  relatively  insensitive  to  the  initial  state. 
However,  we  can  suggest  a  natural  procedure.  Far  from  CPA,  the  aooustic 
azimith  will  change  very  slowly.  Therefore,  we  initialize  the  first  com¬ 
ponent  of  the  state  to  the  first  measurement  received.  The  other  com¬ 
ponents,  which  are  derivatives  of  the  first,  ore  set  to  0.  The  initial  er¬ 
ror  covariance  can  be  chosen  somewhat  more  arbitrarily,  since  error  oovuri- 
anoe  approaches  a  steady  state  value  relatively  quickly,  with  small  values 
of  Q  and  R.  For  the  model  as  determined  above  (R=9 .0,Q=0 .001 ) ,  it  turns 
out  the  covariance  reaches  steady  state  after  approximately  n  measurement 
updates . 

With  the  above  system  model  and  initialization  procedure,  we  can  con¬ 
struct  the  Kalman  filters  necessary  to  implement  the  recursive  probability 
formula  developed  in  the  chapter  4.  This,  in  turn,  gives  us  an  implementa¬ 
tion  for  the  data  association  trees.  We  now  turn  to  the  problem  of  con¬ 
structing  system  models  for  use*  in  the  track  association  filters. 

5.2  -  STATE  SPACE  MODEL  :  TARGET  TRACKING  FILTER 

Developing  a  system  model  that  produces  a  good  filter  to  estimate  the 
full  target  state  from  the  data  tracks  of  the  two  nodes  is  a  much  more  dif¬ 
ficult  task.  Conceptually,  we  can  track  the  target  jtate  in  two  ways.  We 
can  track  the  acoustic  position  and  velocity,  or  we  can  track  the  true  po¬ 
sition  and  velocity.  We  shall  investigate  the  former  method  first. 


X  CmeWs) 


- 


'4 

i 

f 

a 

Consider  Figure  5-5.  Assume  that  the  target  has  u  constant  velocity 
and  constant  heading.  The  node  receives  measurement  6^  ,  at  time  k^, 
corresponding  to  a  true  target  position  at  time  tj.  The  relationship 
between  k1  and  is  given  by  (5.1): 


H1  is  the  range  to  the  acoustic  target  position, 
oeives  the  measurement  «J_,  and  we  have 

c 


Thus, 


+ 


c 


At  time  k2,  the  node  re- 


2  1 

k2-k1  s  t2-t1  +  (5.10) 
Even  if  we  choose  kg-k^  the  sampling  period,  to  be  constant,  the  time 
difference  between  consecutive  acoustic  positions  i3  constantly  changing  in 
a  nonlinear  fashion.  Thus  the  system  equation  is  nonlinear,  even  though 
the  target  trajectory  is  perfectly  linear. 


It  turns  out  that  the  equations  of  the  EKF  for  this  model  are  very 
complicated.  They  turn  out  to  be  very  sensitive  to  linearization.  We  can 
get  an  idea  of  this  by  reasoning  as  follows.  While  the  target  is  approach¬ 
ing  CPA  (during  which,  it  is  hoped,  the  target  will  be  acquired),  the 
acoustic  range  decreases.  According  to  (5.10),  this  means  the  time  differ¬ 
ence  between  acoustic  positions  i3  greater  than  the  time  difference  between 
the  corresponding  measurements.  Now,  since  the  system  equation  describes 
the  time  evolution  of  the  state,  its  nonlinearities  will  affect  the  time 
update  equations  of  the  EKF.  One  would  expect,  then,  that  the  effects  of 
these  nonlinearities  would  be  worse  than  expected  for  a  constant  (known) 
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constant  course),  this  is  not  a  formal  difficulty  as  the  state  variables 
contain  enough  information  to  allow  a  position  to  be  extrapolatod  backwards 
in  time.  Formal  difficulties  arise  when  tho  process  noise  is  assumed  (as 
it  must  be  in  this  case  for  practical  reasons),  but  the  approximation  that 
current  measurements  can  be  derived  from  the  current  state  is  minor  com¬ 
pared  with  other  approximations;  e.g.,  ignoring  altitude,  multipath  signal 
effeots,  etc. 

We  now  turn  to  consideration  of  the  measurement  equation  for  the  sys¬ 
tem  (5.11).  There  are  three  ways  to  incorporate  the  measurements  into  the 
filter.  One  way  would  be  to  combine  the  azimuths  in  some  fashion  to  form 
position  measurements.  This  procedure  is  called  orossfixing ,  and  we  would 
have  a  linear  measurement  equation.  As  discussed  in  the  next  section,  how¬ 
ever,  the  crossfixing  is  not  always  successful,  and  we  could  thus  lose  in¬ 
formation  for  updating  target  state  estimates. 

The  other  two  methods  use  the  given  measurements  from  the  nodes  in¬ 
dependently.  The  first  method  would  identically  initialize  two  filters, 
and  then  run  each  filter  on  a  different  node’s  measurements.  The  estimates 
of  each  filter  could  then  be  combined,  taking  into  account  the  common  ini¬ 
tialization.  The  second  method  would  be  run  one  filter  and  incorporate 
both  measurements.  Theoretically,  both  approaches  should  produce 
equivalent  results.  As  pointed  out  in  the  last  section,  however,  we  do  not 
obtain  much  of  a  decrease  in  the  variance  of  the  position  coordinate  per¬ 
pendicular  to  the  trajectory.  This  implies  that  this  coordinate  i3  rela¬ 
tively  unobservable.  We  wil  thus  obtain  large  errors  in  the  estimation  of 
this  component,  which  will  in  turn  give  poor  linearization.  Without  a  good 
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Since  6=ip-«:,  we  cun  derive  the  following  using  trigonometric  identities: 
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Substituting  into  (5.12)  and  simplifying,  we  get 


t 


~(xvx+yvy)+\!p2o2-(xvy-yvx)2 
l  _____  Z-_ 

C  -V 


(5.17) 


Since  p  and  v  ore  functions  of  x,  y,  v  ,  and  v  ,  (5.17)  gives  t  completely 

x  y 

in  terms  of  the  true  target  state.  We  now  have 
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X-  s  X-tV 
P  X 

yP  =  y-tvy 


and 


-1  XP  -1  x-tvx 

6  s  tdn  t:  =  tdn 


(5.18) 


Thus,  our  measurement  equation  is 


1  x(k)-t(k) vx(k) 

Z<H  *  lun~  y(k)~~tii;)y 00  ’  “lk) 


(5.19) 


where  w  is  the  noise  term  and  t  is  given  by  (5.12). 


Since  this  equation  is  nonlinear,  we  will  have  to  use  an  EKF.  For 
this,  we  need  the  vector  derivative 
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The  equations  for  the  above  partial  derivatives  cun  be  derived  by  straight¬ 
forward  though  tedious,  calculus.  The  results  are 
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where 
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((5.22) 

The  equations  for  the  EKF  arc  of  exactly  the  same  form  as  (4,j),  except  that  here 
the  measurement  matrix  is 
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The  only  thing  left  to  specify  here  is  the  initialization  procedure. 
We  will  initialize  the  position  and  velocity  of  the  filter  at  one  point  in 
time.  This  will  be  done  using  the  smoothed  estimates  of  the  azimuths  and 
their  derivatives  at  each  node.  The  procedure  involves  erossfixing,  which 
is  the  topic  of  the  next  section.  A  poor  initialization  can  cause  problems 
in  two  ways.  If  the  filter  is  initialized  when  the  target  is  far  from  CPA 
of  both  nodes,  the  trajectory  will  be  unaffected  by  the  measurements  for  a 
relatively  long  period  of  time.  The  errors  in  initialization  will  thus 
grow  with  time.  In  addition,  if  the  estimated  target  position  is  much 
closer  to  one  of  the  nodes  than  the  true  target  position,  the  filter  will 
give  greater  weight  to  the  measurements  from  that  node  than  it  should. 
These  effects  cascade,  resulting  in  filter  divergence. 

With  these  points  in  mind,  we  formulate  the  following  restrictions  on 
the  initialization  procedure.  First  of  ail,  we  use  two  azimuths  to  calcu¬ 
late  a  position  only  if  their  difference  i3  greater  than  60°.  This  ensures 
that  we  do  not  initialize  when  the  acoustic  position  is  too  far  away. 
Second,  we  artificially  adjust  the  initial  heading  so  that  it  points  toward 
the  midpoint  of  the  line  segment  joining  the  two  nodes.  This  helps  keep 
the  initial  position  away  from  either  node.  This  procedure,  of  course, 
works  well  only  for  trajectories  that  pass  between  the  nodes.  However, 
this  is  not  much  of  a  limitation  if  we  consider  the  two  nodes  as  part  of  a 
larger  network.  In  such  a  system,  it  is  reasonable  to  assume  thut  the  in¬ 
teresting  targets  will  eventually  puss  between  some  pair  of  nodes.  In  ad¬ 
dition,  u  node  pair  may  receive  an  initial  state  for  a  target  that  is  being 
tracked  by  other  nodes.  This  proocJure,  called  target  hondoff ,  is  impor- 


tant  in  a  system  with  a  large  number  of  nodes.  In  these  systems,  initiali¬ 
zation  is  required  only  if  data  tracks  cannot  be  associated  with  a  priori 
targets.  This  reduces  the" number  of  initializations  a  node  pair  must  per¬ 
form  . 

S.3  -  CROSSFIXING  AND  MON-CLASSICAL  GHOSTS 

As  mentioned  in  the  last  section,  crossfixing  is  the  procedure  for  ob¬ 
taining  position  measurements  from  the  azimuth  measurements  of  the  two 
nodes.  For  time-delayed  measurements,  the  process  is  more  complex  than 
simple  triangulation.  Suppose  that  time  target,  position  is  further  away 
from  node  A  than  node  B.  The  signal  will  thus  reach  B  first.  We  cannot 
obtain  a  position  using  this  measurement  at  B  because  the  corresponding 
measurement  at  A  has  not  yet  arrived.  On  the  other  hand,  if  we  use  the 
current  measurement  of  A,  the  corresponding  measurement  of  B  is  a  past 

measurement.  Thu3,  crossfixing  can  only  produce  acoustic  positions 

corresponding  to  measurements  from  the  node  farthest  away  from  the  target. 
Given  the  current  measurement  from  one  node,  we  can  divide  the  crossfixing 
operation  into  two  steps.  The  first  step  determines  the  measurement  and 
its  time  of  reception  at  the  second  node  that  corresponds  to  the  same 
acoustic  position  33  the  given  measurement .  The  second  step  is  standard 
triangulatiori.  Since  the  first  step  is  the  real  problem  here,  we  will  in¬ 
vestigate  it  hero.  The  derivation  follows  [4], 

Consider  Figure  5-7 .  Suppose  that  we  receive,  at  time  the  azimuth 

<5,,  at  node  B.  We  wish  to  find  the  measurement  «S,  received  by  node  A  at 

b  A 

time  t^  that  corresponds  to  the  same  acoustic  position  :;s  6^.  Let 


6t  a  fc  — t 
B  A 


(5.23) 


ami  define  tp  to  be  the  line  ut  which  the  true  target  position  was  the  same 

as  the  desired  acoustic  pc  ition  p,  Then 

6t  =  (tB-tp)-UA-tp)  (5.24) 

Now,  ty-tp  is  the  travel  time  of  the  signal  between  P  and  B,  while  t,— t  is 

A  P 

the  travel  time  of  the  signal  between  P  cud  A.  From  the  geometry  of  Figure 


5-7,  the  following  is  evident: 


l  _  h _ h 

A  "  sin8A '  "  eosri^ 

l  -  h  _  h 

B  "  sind'  =  cos<5n 

a  b 


(Remember  8^  is  a  negative  quantity  in  this  figure).  Also, 


(5.25) 


(5.26) 


tA“tP  =  c 


(5.28) 


VV  =  r 


From  (5.25)  and  (5.26)  we  get 


(5.29) 


1-1  - « _ 

B  A  sin(8A-iiB) 

Substituting  this  into  (5.27)  and  solving  for  1A,  we  get 

d  cos8g 

1A  '  sin(6A-8B) 

Also , 


(5.30) 


(5.3D 


d  cosd 

1B  ‘  &1  n  ( ) 


(5.32) 


Using  (5.24),  (5.23),  (5.25),  (5.31),  and  (5.32),  we  finally  get 

,  cosd  -  cos8 

6t  a  -  ’  A _ B 

c  '  sinie.  -c..) 

a  b 
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Although  we  have  used  a  restricted  geometry  here,  the  result  is  more  gen¬ 
eral.  The  equation  is  valid,  of  oourse,  only  in  the  range  where  the  meas¬ 
urements  6.  and  <S„  intersect;  namely, 

n  D 


I. 

4  <  6  <  ? 

C  DC 

*  **B 

<  *A 

<  * 

2 

II. 

n 

•  2 

<  *A 

<  21?+^ 

III. 

f  <  i  " 

1  < 

•  2 

*A  < 

^B 

For  6„  =  4,  the  acoustic  position  is  indeterminate. 

13  —cl 

Now,  suppose  that  we  have  the  past  history  of  received  measurements  at 
node  A  (i.c.,  the  azimuth  V3  time  curve  for  A)  us  in  Figure  5-b.  Using 

(5.32)  we  c>,n  convert  the  single  observation  dj.  ax.  node  B  into  a  curve  of 

possible  6^  vs  6t  and  plot  it  on  the  sane  scale- (Figure  5-y).  This  curve 

is  called  the  reflected  observation  curve .  The  intersection  of  these  two 

curves  gives  the  desired  measurement  and  its  reception  time  t^. 

It  turns  out  that  multiple  intersections  of  the  two  curves  is  possi¬ 
ble.  This  occurs  when  the  target  trajectory  passes  between  the  two  nodes. 
Figure  5-10  shows  an  example.  At  time  instants  one  second  apart  we  have 
taken  the  noiseless  measurements  at  13,  and  obtain  .ill  possible  acoustic  po¬ 
sitions  assuming  a  full  past  history  of  noiseless  measurements  at  A.  The 
target  velocity  whs  150  meters  per  second  at  a  heading  of  135  degrees  from 
north,  and  the  distance  between  A  and  B  is  5000  meters.  In  the  figure,  the 
real  acoustic  trajectory  is  obvious.  The-  false  positions  form  a  track  that 
flies  from  the  perpendicular  bisector  of  the  line  segment  joining  A  and  B 

—  J  *“ 


I 
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to  node  A.  This  behavior  is  exhibited  by  all  such  false  trucks,  and  they 
occur  only  on  the  incoming  path.  We  call  those  tracks  non-classical 
ghosts .  to  distinguish  then  from  the  ghosts  of  Chapter  4. 

A  heurestic  procedure  can  be  developed  which  eliminates,  for  noiseless 
measurements,  most  traces  of  ghost  tracks.  Looking  again  at  Figure  5-10, 
we  see  that  the  ghost  track  is  concave  downward.  (The  track  begins  at  the 
biseotor  and  heads  toward  A).  Since  the  target  heading  is  135  degrees,  the 
acoustic  azimuth  increases  with  time.  However,  succeeding  positions  along 
the  ghost  track  requires  the  acoustic  azimuth  to  decrease  with  time.  Thus, 
the  measurements  of  A  that  determine  the  ghost  track  forms  a  sequence  that 
goes  backwards  in  time.  At  each  point  in  time,  then,  after  crossfixing  has 
been  completed,  we  eliminate  the  past  history  at  A  that  occurs  prior  to  the 
earliest  crossfix,  then  the  rest  of  the  ghost  track  will  be  eliminated.  If 
it  is  not,  the  ghost  track  remain  .  Figure  5-11  shows  the  results  for  this 
example.  As  can  be  seen,  only  part  of  the  ghost  track  is  eliminated. 

We  can  eliminate  the  rest  of  the  ghost  by  observing  that  this  portion 
of  the  ghost  track  must  correspond  to  the  later  crossfix,  else  it  would 
have  been  eliminated.  We  thus  keep  only  the  earliest  crossfixes,  which 
correspond  to  tne  real  acoustic  positions.  The  result,  with  a  few  minor 
adjustments  that  need  not  concern  us  here,  are  shown  in  Figure  5-12,  Mote 
that  the  procedure  makes  one  expected  error  —  the  first  point  at  which  the 
ghost  is  the  earliest  crossfix. 

In  the  noisy  case,  deghosting  is  somewhat  more  difficult,  primarily 
because  crossfixing  is  not  alwuy3  successful.  An  important  thing  to  notice 
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in  Figure  5-12  is  the  gap  that  occurs  in  the  region  where  the  ghost  traok 
and  the  target  track  intersect.  In  this  region,  the  reflected  observation 
curve  is  almost  tangential  to  the  azimuth  time  curve;  the  intersection 
points  are  close  together.  This  is  shown  in  Figure  5-13  .  In  cases  like 
these,  it  is  very  difficult  to  detect  the  intersection  points.  This  is  why 
the  gap  appears  in  Fig.  5.12;  the  curves  were  almost  tangential  and  the 
crossfix  failed.  The  situation  is  even  worse  when  the  measurements  are 
noisy.  The  gaps  can  be  many  time-sampling  periods  long.  In  addition  to 
this,  the  noisy  data  can  cause  random  crosfixirig  failures  anywhere  along 
the  path.  Because  of  thia  loss  of  information,  we  choose  not  to  use  posi¬ 
tion  measurements  as  updates  to  our  filters.  We  will,  however,  U3e 
crossfixing  in  initialization. 

5.4  -  SUMMARY 

In  this  chapter,  we  have  developed  stats  variable  descriptions  neces¬ 
sary  for  the  construction  of  Kalman  filters  for  both  the  data  association 
trees  and  track  association  trees >  We  also  discussed  the  operation  of 
crossfixing  and  non-classical  ghosts.  In  chapter  6,  we  shall  present  some 
results  for  this  implementation  of  our  tracking  algorithm. 
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CHAPTER  6 

RESULTS 

In  this  chapter,  we  present  demonstrations  of  the  filters  developed  in 
Chapter  5  and  their  application  to  the  tracking  algorithm.  It  should  be 
pointed  out  that  these  examples  only  highlight  basic  characteristics.  In 
order  to  evaluate  the  strength  of  the  filters,  Monte  Carlo  simulations 
should  be  conducted  over  a  wide  range  of  target  scenarios,  and  performance 
statistics  must  be  defined  and  computed  to  facilitate  the  evaluation.  Un¬ 
fortunately,  due  to  time  constraints,  we  were  unable  to  give  a  really 
thorough  performance  evaluation.  We  will,  however,  present  some  results 
that  illustrate  the  characteristics  of  the  filters  and  the  algorithm. 

All  of  the  demonstrations  in  this  chapter  were  done  by  computer  simu¬ 
lation.  The  simulations,  us  well  as  the  algorithm  itself,  was  written  in 
the  C  language  and  run  on  a  PDP-11  computer. 

In  the  first  example,  a  Monte  Carlo  simulation  was  run  on  the  data  as¬ 
sociation  filter  to  provide  some  information  on  its  performance.  The 
parameters  of  the  example  trajectory  are  those  of  Target  1  in  Table  6-1. 
The  node  A  at  which  the  data  association  is  carried  out  is  positioned  at 
the  origin.  A  set  of  noiseless  acoustic  azimuths  was  generated  for  A,  as¬ 
suming  a  sampling  period  of  one  second.  Only  azimuths  that  corresponded  to 
angles  around  CPA  that  were  between  -70°  and  70°  were  retained.  From  this 
azimuthal  set,  lj0  sets  of  noisy  data  were  generated  by  adding  white  Guas- 
sian  noise  with  a  standard  deviation  of  j°.  The  data  association  filter 


was  run  separately  on  each  set  of  data.  For  each  point  in  time,  an  average 
mean  square  error  was  calculated  over  the  '30  data  sets.  The  results  are 
shown  in  Table  6-2.  Clearly,  the  filter  performance  is  best  far  from  CPA, 
when  the  azimuth-time  curve  is  practically  linear.  As  the  acoustic  data 
approaches  the  CPA  angle,  however,  the  azimuth  derivatives  are  in  flux,  and 
performance  degrades.  Still,  the  overall  average  squared  error  was  reason¬ 
able,  about  half  the  variance  of  the  measurement  noise. 

The  next  example  is  an  excellent  demonstration  of  the  capabilities  of 
the  data  association  process  in  a  two-target  environment.  The  trajectories 
are  those  of  Target  1  and  Target  2  in  Table  6-1.  Again,  the  reference  node 
A  is  at  the  origin.  Figure  6-1  shows  the  exact  acoustical  data  generated 
by  these  trajectories.  Note  that  the  curves  intersect  twice,  and  that  the 
angles  of  Intersection  are  somewhat  sm-Jll,  Tnsse  regions  are  potential 
trouble  spots  for  the  data  association  process  in  the  presence  of  noise. 
The  node  must  be  able  to  distinguish  tue  noisy  data  so  that  they  can  be 
tracked.  Otherwise,  the  target  tracking  process  will  not  be  good.  Figure 
6-2  shows  the  noisy  data  after  the  data  association,  while  Figure  6-3  plots 
the  filtered  estimates.  Although  a  few  measurements  around  the  curve  in¬ 
tersections  points  are  mixed  up,  the  process  manages  to  lock  onto  and  track 
separately  the  two  3ets  of  data . 

The  next  example  demonstrates  the  target  tracking  filter.  A  second 
node  B  is  placed  at  the  point  ('3000,  0).  Noisy  sets  of  data  are  generated 
for  both  A  and  B  from  the  trajectory  of  Target  1.  The  resulting  target 
track,  combining  the  two  data  tracks,  is  shown  in  Figure  6-4,  Because  of 
the  restrictions  imposed  in  Section  *3.2,  the  target  was  well  along  the 
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trajectory  before  the  track  was  initialized.  The  plot  is  a  bit  misleading, 
in  that  the  points  of  the  track  do  not  lie  in  sequential  order  following 
the  true  trajectory.  This  because  the  estimate  of  the  current  target  state 
is  updated  by  current  measurements,  which  correspond  to  earlier  states  of 
the  target.  There  always  exists  an  uncertainty  in  the  time  of  occurence  of 
the  earlier  states;  hence,  one  would  expect  consecutive  estimates  would  not 
necessarily  follow  a  sequentially  follow  the  true  trajectory. 

Our  last  example  demonstrates  track  association  with  the  data  tracks 
of  the  second  example.  Figure  5-5  shows  the  results.  The  ghost  tracks 
(i.e.,  the  incorrect  truck  associations)  turned  out  to  be  not  much  of  a 
problem.  One  ghost  did  not  even  initialize  a  track.  The  other  ghost  in 
due  course  surpassed  the  speed  of  sound,  and  the  algorithm  automatically 
rejected  it. 

The  peculiar  behavior  of  the  track  of  Target  2  is  easily  explained  in 
terms  of  the  initialization  procedure  outiineJ  in  chapter  5.  The  initial 
trajectory  was  set  to  head  towards  the  point  (2500,0).  The  track  more  or 
less  held  this  direction  until  the  target  approached  CPA.  Then  the  meas¬ 
urements  from  both  nodes  began  to  have  effect,  and  the  track  was  subse¬ 
quently  pulled  back  around  the  actual  trajectory. 

As  we  mentioned  earlier,  the  examples  only  highlight  the  major  charac¬ 
teristics  of  the  tracking  algorithm  in  the  two-node  system.  We  did  provide 
some  statistics  on  the  performance  of  the  data  association  filter,  pointing 
out  its  weaknesses,  and  demonstrated  its  success  in  resolving  long  term  am¬ 
biguities.  We  also  demonstrated  the  full  target  state  filter  and  the  track 


association  process.  However,  much  work  still  remains  to  be  done  in  order 
to  fully  test  the  tracking  algorithm.  Heasureiaents  of  performance  of  the 
full  state  tracker  need  to  be  evaluated.  Also,  the  performance  of  the  data 
association  filter  as  a  function  of  the  rate  of  false  alarms  should  be  in¬ 
vestigated  tc  provide  a  measure  of  robustness. 

The  application  of  these  ideas  to  larger  distributed  networks  is  also 
a  further  area  of  .-a;  -oh  that  must  be  investigated.  This  would  probably 
require  a  theoretical  structue  that  combined  various  local  target  tracks 
into  global  tracks  suitable  for  the  users  of  the  tracking  system.  The  ef¬ 
fect  of  communications,  problems  with  maneuvering  targets  and  target  han- 
doff  procedures  are  other  areas  that  will  require  careful  research.  Final¬ 
ly,  we  mention  the  importance  in  testing  the  algorithm  in  a  real  world  sys¬ 
tem  to  investigate  its  real  utility. 

What  we  have  attempted  here  is  to  show  that  our  tracking  algorithm 
does  perform  reasonably  well  in  tracking  multiple  targets  and  resolving  am¬ 
biguities  in  data  associations.  It  is  hoped  that  it  will  be  a  useful  tool 
for  future  investigation  into  distributed  processing  tracking  systems. 
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